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We discuss theoretically and numerically the intramolecular form factor F(q) in dense 
polymer systems. Following Flory's ideality hypothesis, chains in the melt adopt Gaussian 
configurations and their form factor is supposed to be given by Debye's formula. At striking 
variance to this, we obtain noticeable (up to 20%) non-monotonic deviations which can be 
traced back to the incompressibility of dense polymer solutions beyond a local scale. The 
Kratky plot (q 2 F(q) vs. wavevector q) does not exhibit the plateau expected for Gaussian 
chains in the intermediate g-range. One rather finds a significant decrease according to the 
correction S(F~ 1 (q)) = q 3 /32p that only depends on the concentration p of the solution, 
but neither on the persistence length or the interaction strength. The non-analyticity of the 
above q 3 correction is linked to the existence of long-range correlations for collective density 
fluctuations that survive screening. Finite-chain size effects are found to decay with chain 
length N as 1/y/N. 
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I. INTRODUCTION 

Following Flory's ideality hypothesis JJ|, one expects a macromolecule of size N, in a melt 
(disordered polymeric dense phase) to follow Gaussian statistics . The official justification 

of this mean-field result is that density fluctuations are small, hence negligible. 

Early Small Angle Neutron Scattering experiments [(], Q] have been set up to check this central 
conjecture of polymer physics. The standard technique measures the scattering function S(q) (q 
being the wave-vector) of a mixture of deuterated (fraction f) and hydrogenated (fraction 1 — /) 
otherwise identical polymers. The results are rationalized [f], Q] via the formula 

S(q) cc f(l - f)F(q) (1) 

to extract the form factor (single chain scattering function) F(q). To reveal the asymptotic behavior 
of the form factor for a Gaussian chain F^(q) ~ 12/q 2 b 2 , one usually plots q 2 F(q) versus q 
(called "Kratky plot"). The aim would be to show the existence of the "Kratky plateau" in the 
intermediate range of wave- vectors < q < ^ ("Kratky regime") where R g is the radius of 
gyration of the macromolecule and b is the (effective) statistical segment length [3|. In contrast 
to the low-g "Guinier regime" (q J$ tt-)> clean scattering measurements can be performed in the 

Q g 

Kratky regime [9] suggesting the measurement of b from the height of the Kratky plateau. 

Surprisingly, this plateau appears to be experimentally elusive as already pointed out by Benoit 
[(J: "Clearly, Kratky plots have to be interpreted with care". For typical experiments, the available 
g-range is 5-10~ 3 /A < q < 0.6/A. Kratky plots are quickly increasing at high q because of the 

rod- like effect starting at q ~ 7-, when the beam is scanning scales comparable with the persistence 
length l p (this regime is in fact used to assess l p ). Sometimes these curves can also quickly decrease, 
this is usually attributed to the fact that the chain cannot be considered as infinitely thin. The 
finite cross section of the chain tends to switch off the signal [f], Q] as exp (— q ^ c , where R c is the 



radius of gyration of the cross section (in CSo, a good solvent for dilute polystyrene, R 2 = 9.5A 2 
[10l]). Rawiso et al. [6|, [Z( have also shown that sometimes these two effects can compensate, for 
instance in a blend of hydrogenated and fully deuterated high molecular weight polystyrene, letting 
appear a pseudo Kratky plateau, extending outside the intermediate regime to higher g-values. In 
fact, taking for instance the case of polystyrene with l p ~ 10A, all these parasitic effects never 
really allowed neutron scattering experiments to confirm the Flory's hypothesis. Up to now no 
scattering experiment has been performed on a sample allowing for a test over a wide enough range 



of q 



11 



121 ] and there exists no clear experimental evidence of the Kratky plateau expected for 



3 



Gaussian chains. 

As we will show in this paper, there are fundamental reasons why this plateau may actually 



never be observed, even for samples containing very long and flexible polymers. Recent 



Yjlong- 



range correlations, induced by fluctuations, have been theoretically derived |13l. 1 141. Il5l. Il6l. 1 171. 116 1 
and numerically tested for two JjJ] and three-dimensional [itJ dense polymer systems. (Similar 
deviations from Flory's ideality hypothesis have been also reported in various recent numerical 

221 . |23j.) The conceptually simpler part of these 



studies on polymer melts 120, 211 and networks 

Q 

effects is related to the correlation hole [2j and happens to dominate the non-Gaussian deviations 
to the form factor described here 2^ |. 

This is derived in the following Section [XTJ. There we first recapitulate the general perturbation 
approach (Sec. Ill A|) . discuss then the intramolecular correlations in Flory size-distributed poly- 
mers (Sec. Ill Bp . We obtain the form factor of monodisperse polymer melts by inverse Laplace 
transformation of the polydisperse case (Sec. Ill Cj) . In Section ITTT1 these analytical predictions are 
illustrated numerically by Monte Carlo simulation of the three dimensional bond-fluctuation model 
26]. We compare melts contain ing only monodisperse chains [jj], [it]] with systems of (linear) 



25, 



equilibrium polymers (EP) 



27 



28, 



29( where the self-assembled chains (no closed loops being al- 
lowed) have an annealed size-distribution of (essentially) Flory type jsfj. Excellent parameter free 
agreement between numerical data and theory is demonstrated, especially for long EP (Sec. IIIIB]) . 
Finite-chain size effects are addressed as well. In the final Section HVl the experimental situation is 
reconsidered in the light of our analytical and computational results. There we show numerically 
that eq[U remains an accurate method for determining the form factor which should allow to detect 
the long-range correlations experimentally. 
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II. ANALYTICAL RESULTS 



A. The Mean-Field Approach 



It is well accepted that at the mean-field level [3] , the excluded volume interaction is entropically 
screened in dense polymeric melts. Long ago Edwards and de Gennes 



32, 



331 ] developed 



a self-consistent mean-field method to derive a screened mean (molecular) field: this theory is an 



adaptation of the Random Phase Approximation (RPA) [34J] to polymeric melts and solutions. The 
famous result of this approximation gives the response function S(q) as a function of F^°\q), the 
scattering function of a Gaussian (phantom) chain via the relation: 

1 _ 1 

where p is the mean concentration of monomers in the system and v is the bare excluded volume 
(proportional to the inverse of the compressibility of the system). Please note that F^(q) has 
to be properly averaged over the relevant size-distribution of the chains [33]. To get the effective 
interaction potential between monomers, we label a few chains. The interactions between labeled 
monomers are screened by the background of unlabeled monomers. Linear response gives the 
effective q dependent excluded volume: 

_J_ = i + pF(°)(g). (3) 

Let us from now on consider a dense system of long chains with exponentially decaying number 
density p^ = pp 2 e~ flN for polymer chains of length N with u = 1/ (iV) being the chemical potential. 



28]. Hence, eq[3]yields (using eqs[6 



This so-called Flory distribution is relevant to EP systems 
indicated below) 

fi g >£ 1 + q z £ z a z 

Here a is the characteristic length of the monomer (a 2 = b 2 /6) and f = yjl^ is the mean-field 
correlation length. When chains are infinitely long (p — > 0) we recover the classical result by 
Edwards [3] ignoring finite-size effects. If we further restrict ourselves to length scales larger than £ 
(g£ <C 1) eq H] simplifies to v e ff(q) « q 2 a 2 /2p which does not depend on the bare excluded volume v 
and corresponds to the incompressible melt limit. For very large scales (qR g <C 1) one obtains the 
contact interaction associated to the volume v* = v e g(q — ► 0), such that v* /v = p£, 2 /a 2 ~ £, 2 /Rg 
(far weaker than the initial one given in the direct space by v(r) = v5(r)). The interaction v* is 
relevant to the swelling of a long chain immersed in the polydisperse bath Q]. 
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The screened excluded volume interaction eqH] taken at scale R g is weak and decays with chain 
length as 1/(N). The associated perturbation parameter u in d-dimensional space depends on 
chain length as u ~ v*N 2 /Rf ~ (N) 1 ~ d ^ 2 and the screened excluded volume potential is, hence, 



a. 



perturbative in three dimensions 

Let us define G(q, s) = (exp(— iq • r s )), the Fourier transform of the two-point intramolecular 
correlation function, with s = the number of monomers between the two positions separated 

by r s . One can perturb the two-point Gaussian correlation function 

(q, s) = <exp(-iq • rj) = exp (sq 2 a 2 ) (5) 

with the molecular field eq[H In this type of calculations, there are only three non-zero contributions 



35, 



3a]. They are illustrated by the diagrams given in Figure Ufa). Knowing this correlation, it is 



possible to derive many single-molecule properties 



B. Intramolecular correlations for Flory size-distributed polymers 



The intramolecular correlation function is investigated through its Fourier transform, the form 
factor. As already mentioned, we consider Flory size-distributed polymer systems [30]. In this 
case, one can define the form factor as: 

oo N N 

F(q) =J2^ e*P(-/^0 E E< ex PHq ' r u)> ( 6 ) 

Af=0 i=l j=l 

with q the scattering vector. If the chains followed Gaussian statistics (as suggested by the Flory's 
hypothesis), one should find using eq[5j 



1 



(7) 



H 1 + q 2 a 2 /n 

The ideal chain form factor for Flory size-distributed polymers is represented in the Figures El [3] 
and H] where it is compared to our computational results on EP discussed below in Section ITO1 In 



the smal 



-q regime, qR g <C 1 and for a polydisperse system, we can measure R g via the Guinier 



relation \§, la, 



33J: 



Jn(0)2 

1 \z 



(8) 



Please note the Z-averaging 



331 ] in the definition of R g) z 

Z N N 2 PN R 2 JN) 



U R,2 



EnN 2 Pn 



0) 
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FIG. 1: (a) Interaction diagrams and (b) integration contour and singularities, (a) There exist three 
non-zero contributions to first-order perturbation, the first involving two points inside the segment of size 
s = \i — j\ monomers, the second involving one point inside and one outside the segment and the third 
involving one point on cither side of the segment. Momentum q flows from one correlated point to the other. 
Integrals are preformed over the momentum k. Dotted lines denote the effective interactions v e ff(k), bold 
lines the propagators, (b) A possible contour in the complex plane. The dashed line depicts the logarithmic 



branch-cut, which is an arc of the circle of radius r = yji + q 2 a 2 between q\ = re 1 ^ a > and <?2 = re la , 
with a = arctan(^/ n/q 2 a 2 ). The other singularities are single poles at k = k = iyj n/a 2 and at 



k = i\/ q 2 + [i/a 2 . 

Rg(N) being the radius of gyration in the monodisperse case. (Since only Z-averaged length 
scales are considered below the index Z is dropped from now on.) As the Flory distribution gives 
(N p ) = T(p+ one has Rf )2 = 3a 2 //"- In the Kratky regime between coil and monomer size 

one recovers the classical result for infinite chains 

2 



F(°)(g) 



(10) 



(indicated by the dashed line in Figure E|) which expresses the fractal dimension of the Gaussian 
coil. 

Perturbing the Gaussian correlation function with the screened potential u e ff( r ) = 
/ - yefT(<?K q ' r d 3 q/(27r) 3 , one has to evaluate (SA) = - (UA) + (U} (A) with U = \ j Vgsfa-rj) 
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FIG. 2: The intramolecular form factor F(q) of polymer chains is an important property since it allows 
to make a connection between theory and simulation on the one hand and experiments of real systems on 
the other hand. The ideal chain form factor for Flory size-distributed polymers, eqlTJ is indicated by the 
solid line. In the Kratky regime between the total chain and monomer sizes the form factor expresses the 
fractal dimension of the Gaussian coil, eq [TO] (dashed line) . Experimentally, this is the most important 
regime since it is, for instance, not affected by the (a priori unknown) polydispersity. The form factors of 
equilibrium polymers (EP) of various scission energies E obtained numerically are indicated. The form factor 
of polydisperse polymer systems is obtained by computing for each chain Q]\ sin(q • rj) 2 + (^\ cos(q • r 4 )) 2 , 
summing over all chains (irrespective of their length) and dividing by the total number of particles. The 
computational data reveal an additional regime at wavevectors corresponding to the monomer structure 
("Bragg regime") which is not treated by our theory. All data have been obtained for a number density 
p = 0.5/8 of the three dimensional bond-fluctuation model (BFM). 

for the observable A of interest Q]. In our case, A corresponds the form factor, eq[6] (without aver- 
age) . The calculation in reciprocal space is schematically illustrated by the three diagrams given in 
Fig. DJa) where bold lines represent propagators and dotted lines the effective interactions v e s(k). 
For the first perturbation contribution, {UF) , the propagator carrying a momentum p (with 
p = qorp = q — kas indicated in the diagrams) corresponds to a factor l/(p 2 a 2 + /u). For the 
second contribution, {U) {F) , the same rules apply but q has to be replaced by zero and (then) 
/j, by fi + q 2 a 2 . The momentum k is integrated out. Each diagram gives a converging contribu- 
tion. Summing up the contributions of the diagrams (the central one has to be counted twice) 
and performing the angular integral (over the angle between k and q) we arrive at the following 
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integral: 

5F(q)=F(q)-F^(q) = 

v f ,, k*e (2^ + (q 2 + k 2 )a 2 ) 2 

4Tr 2 a 2 (fi + q 2 a 2 yJ R l + k 2 e fM + k 2 a 2 I 11 ) 



2 1 , ///+ (q + k) 2 a 2 

in 1 



v /x + (q 2 + A; 2 )a 2 2kqa 2 \fi + (q — k) 2 a 2 
The contributions to this integral come from two poles, one at k = , this high-A; contribution 



renormalizes the statistical segment, and one at k = iy q 2 + and from the logarithmic branch 
cut. The integration contour and the singularities are illustrated in Figure [2(b). Absorbing the 



high-/c pole contributions in the renormalized statistical segment b* 



jj], and using b* instead of b 

(°) - q™2 



in the definition of F^(q), one finds (in the limit q£ < 1) as a function of Q = qRg = 3qa 2 / \x 

SF(Q)/c = ^ ( - 1 )- i arctan (-9-) + ± 2 + & , (12) 

Q 2 ^(1 + ^1)3/2 J Q \2V3J V3 (1 + ^)2 

gR (p) 

We have introduced here the factor c = — ^ to write the deviation in a form which should scale 
with respect to chain length. The statistical segment length b* , we have introduced, is given by 

V* = 6 2 (l + i^) . (13) 



37|, 



38]. 



which is consistent with the result obtained by Edwards and Muthukumar 

As may be seen from the plot of eq[l2]in Figure[5]the deviation from ideality is positive for small 
wave- vectors (with a pronounced maximum at Q ~ 1). It becomes negative when the internal coil 
structure is probed (Q> 2). Asymptotically, eq [T21 gives 

1 1 

SF(q) ~ -^-7? 2 4 0)2 , (14) 
(dashed line in Figure [5]) which highlights the average swelling factor of the molecule in the melt: 

(comparable to a logarithmic term in the two-dimensional case 13j, [39( ) . This is a sign of swelling, 
because R g = (N) 1 / 2 f(y/{N)), with /, an increasing function, showing that the apparent swelling 
exponent v for finite (N) is slightly larger than 1/2. It is interesting to compare it with the (also 
Z-averaged) end-to-end distance, easily available because it involves only the top-left diagram of 
Figure [TJ 
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FIG. 3: Successful scaling of the form factors of EP obtained from our BFM simulations for various scission 
energies E as indicated: F(Q)/F(0) is plotted vs. the reduced wavevector Q = R g q. Note that both scales 
F(Q -4 0) = (N 2 ) I (N) and the (Z -averaged) gyration radius R g have been directly measured for each 
sample and are indicated in Table HI Obviously, the scaling breaks down due to local physics for large 
wavevectors (Bragg regime). The line represents the prediction for ideal Flory-distributed polymers, eqEl 
with parameters chosen in agreement with the Guinier limit, eq[8] Importantly, in the intermediate Kratky 
regime small, albeit systematic, deviations are visible which will be further investigated below. 



(The diagram must be twice differentiated with respect to q.) The naively denned size dependent 
effective statistical segment of an iV-chain (from B% = b^ s N) therefore is: 

const 



bin 



1 



pb* 3 VN, 

The size-dependences in eqs 1151161 follow the same scaling, but the numerical factors are different. 



(17) 



Although internal segments carry a smaller correction [l4 |. the size-dependent contact potential 
p,/2p in eq. [4] counterbalances this effect and makes the correction to eq[l6]a little smaller than 
the one in eq [T5l 



The asymptotic behavior of eq[12]in the Kratky regime gives 

12 / 3^/6 ^ 3q 



SF(q) 



u*2 ' 



(18) 



<zRg>i q 2 b* 2 yTTpb* 6V ' 8pb* 

which is represented in Figure [5] by the dashed-dotted line. The first term in this equation is a 
size-dependent shift of the Kratky plateau, and the second one, independent of the size, makes 
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FIG. 4: Kratky representation of the intramolecular form factor F(q)q 2 vs. wavevector q for monodis- 
perse (crosses) and equilibrium polymers. The non-monotonous behavior predicted by the theory is clearly 
demonstrated. The ideal chain form factor, eq [7] (thin line), overpredicts the dip of the form factor at 
q « 0.7 by about 20%. The bold line indicates the prediction for infinite chains, eq[23l which should hold 
for both system classes for infinitely long chains. For this reason we have chosen the largest chains currently 
available. 



the essential difference with the Flory prediction (bold line in Figure Hence, the corrections 
induced by the screened potential are non-monotonic (Figure [5j) . 

Eq [11] is not restricted tojjl; < 1 , it is applicable over the entire g-range in the case of weakly 
fluctuating dense polymers 



401 ] (mean-field excluded volume regime) as may be simulated with soft 
monomers allowed to overlap with some small penalty. In the case of strong excluded volume and 
less dense solutions (critical semidilute regime, not explicitly considered here) the results are valid 
at scales larger than £ provided the statistical segment is properly renormalized. Quantitatively, 



the Ginzburg parameter measuring the importance of density fluctuations reads G\ 



;/( P 6 6 ). 



For persistent chains l p > t, t being the thickness of the chain, density fluctuations are negligible 



provided G z ~ 1/ '{pl^/t ) -C 1, with p = pt , the monomer volume fraction. The above makes sense 
if £S> l p , which requires p < t/l p . This criterion also indicates the isotropic/nematic transition 

3 



41]. In summary, mean- field applies provided < p < j-. 
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FIG. 5: Correction SF(Q) = F(Q) - F (0 \Q) to the form factor as a function of Q = qR { g } (using ^ = 
3a* 2 / jj) as predicted by eq[T2]for Flory-distributed polymers. The deviation is positive for small wavevectors 
(Guinier regime) and becomes negative at about Q i=s 3. The scaling factor c = 9R^ j-nph*^ should allow a 
data collapse irrespective of the mean chain size — provided that the chains are sufficiently long to suppress 
additional physics in the Bragg regime. Also included are eq Q3] and [TS] for the asymptotic behavior in 
the Guinier and the Kratky regimes respectively. The data from our BFM simulations (given here for two 
scission energies where high precision data are available) agree quantitatively (especially for small Q) with 
eq HU (Note that the chains are too short to allow a better fit for large Q.) The reference form factor 
F^\Q) has been computed from eq [7] supposing a perfect Flory distribution. 



C. Monodisperse polymer melts in three dimensions 



It is possible to relate the form factor of the polydisperse system (Flory distribution) to the 
form factor F^(q) of a monodisperse system. Following eq[6j 



EV \ Y,N=QPNNF N (q) \^ 2 ( AT\ AT TP t \ 



(19) 



N=0 



the deviations of the form factors of monodisperse and polydisperse systems are related by the in- 
verse Laplace transformation 5Fj\f(q) = jjC^ 1 I -^5F(q) I , C being the Laplace transform operator. 



Using our result eq[T2]ibr polydisperse systems this yields: 



SF N (Q)/c 



+ 



V^Q 2 



V^Q 3 



+ 



V^Q 2 \V^Q V^Q 3 
Q 



D(Q) 



2 faQ 

3 V 2 



ft- 



da, 



(20) 



12 
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+ 



Eq 20 

Eq 21 : N=2048 
Eq 22: N=2048 
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N=256 
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N=1024 
N=2048 
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FIG. 6: Deviation to the ideal form factor in a monodisperse system, eq[20] The shape is similar to that 
in the polydisperse case (Figure O since polydispersity merely affects the Guinier regime. Also indicated 
are the asymptotic behavior in the Kratky regime, eq I21i and the prediction for infinite chains, eq I22l The 
simulation data scales roughly, especially in the Guinier limit (hump) , while for larger q deviations are visible 
which may be attibuted to local physics. The agreement with eqJ^Olis only qualitative here since a much 
more proncounced non-monotonic behavior is seen in the simulation. This is due to use of the polydisperse 
chain perturbation potential, eq|4] — an approximation which must become insufficient for low wavevectors 
where the coil size matters. 



where Q = qR^ N with R 



M 

L g,N 



Na, the radius of gyration in the ideal monodisperse case and 



D(x) is the Dawson function, whose definition is D(x) = e~ x2 f^e^dt 42]. Eq I20I shown in 
Figure [6] (bold line), is accurate up to the finite size corrections to the interaction potential as 
these have been calculated for the Flory distribution, eq[U However, on small length scales, this 
influence is weakened, and in this limit, it gives: 

6^6 



F(q) 



12 



<?R g >i q 2 b 



1 + 



1 



3q 



ir 3 / 2 pb* 3 ^/N 



* 2 



(21) 



Eq[2T]is represented in Figure [6] by the dashed-dotted line [47|. Taking (N) = eq [181 deviates 
from eq [21] only by the numerical coefficient in front of l/y/N. The difference is ~ 10%. 
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D. Infinite chain limit and scaling arguments 

It is worthwhile to discuss the infinite chain limit, p, — ► 0, that puts forward most clearly the 
essential differences with an ideal chain. (The differences between monodisperse and polydisperse 
systems are inessential in this limit.) We can write the form factor of an infinite chain (i? g —* oo) 
at scales larger than £ as: 

F (Q) = -^17 ( 22 ) 



q 2 b* 2 (i + 



Following standard notations 



19, 



43 



44 ], we may rewrite eq[22]in the form: 
1 q 2 b* 2 1 q 3 



(23) 



F(q) 12 32 p ' 

See the Figures H] (bold line), [6] (dashed line) and [9] (bold line) for representations of this important 
limiting behavior. The correction term obtained in the one-loop approximation (eq 123ft depends 
neither on the excluded volume parameter v nor on the statistical segment. Hence, it is expected to 
hold even in the strongly fluctuating semidilute regime and it is of interest to compare our results 
with the recent renormalization group calculations of L. Schafer [43J]. There the skeleton diagrams 



for the renormalization of interaction and statistical segment have also been performed within the 
one loop approximation. From the above it is expected that both results are identical. After careful 
insertions (eq 18.23, p. 389 of Ref. [43j) a q 3 /p correction can be extracted with the universal 
amplitude 0.03124... The fact that this numerical amplitude is so close to our 1/32 = 0.03125 
comforts both our and Schafer's result. 

Performing an inverse Fourier transform of the form factor, eq 1221 gives not only the Coulomb- 
like term from the singularity at the origin, but also another long-range contribution arising from 
the branch cut: 

12 1 _ 9 1 1 

The correction is never dominant in real space. But both contributions are different in nature. In 
the collective structure factor [l(| S~ 1 {q) ~ v + c<iq 2 + c^q 3 , the leading singularity of eq |2"21 
is shifted away from the origin and the corresponding contribution is screened on the lengthscale 
£ in real space. The branch cut (from the q 3 term) still contributes a power law, namely an anti- 
correlation term decreasing as r~ 6 , that has been identified as a fluctuation- induced Anti-Casimir 



effect [15I. [l6l| (or as the n — 1 = —1 Goldstone mode in the polymer-magnet analogy 45]). The av- 



erage number of particles from the same molecule in a sphere of radius R, M{R) = J Gi n t Ta (r)d 3 r 
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is decreased (compared to a Gaussian coil) because of the sign of the correction. Nevertheless, 
the differential (apparent) Hausdorff dimension [46J du as defined by dn = d log M(R) /d log R is 
increased. 

The fluctuation corrections presented in eqs 1221231241 can be interpreted with the following 
argument, involving G(r,s), the correlation function of two points separated by s monomers, 
G(r,s) = f G(q, s)e' q ' r . For large s, the difference SG(r,s) = G(r,s) — &°^(r,s) is discussed 
below in the limit of large and small geometrical separation r between monomers. Here G^(r,s) = 
const exp ( — | J , where the renormalized b* is used instead of b. 



2 sb*\ 

A highly stretched s-fragment (r 3> b*yfs) can be viewed as a string of Pincus blobs, each blob 
of g ~ s 2 b 2 jr 2 units. Different blobs do not overlap in this limit, therefore the effective statistical 
segment of the blob, b e s comes as a result of interactions of units inside the blob (see eq[T7|) 

*- i -( I -^f) 

where const is a universal numerical constant. The elastic energy of stretching the s-segment is 
therefore 

3(r/n b ) 2 3 r 2 3 r 2 , 
W e i(r) ~ n b - — -= — = -— =- ~ -—j 1 + const 



'2 gb 2 s 2sb 2 cS 2 s b* 2 V pb* A S/ 
where n& = s/g is the number of blobs in the s-segment. Thus, when r 3> y^sfc*: 

G(r, s) = const e ~ w ^ ~ G (0) (r, s) exp ^-const ^ (26) 

The faster decay of G as compared to d?( ) leads to higher scattering (positive 5F(q)) at small q. 

On smaller length scales, r <C yfsb* , it is convenient to consider the s-fragment as a chain of blobs 
of size r, with g ~ r 2 /b* 2 . This is sketched in Figure [71 The correction SG here is essentially due 
to the direct interaction of the overlapping blobs or radius r around the two correlated monomers. 
The number of binary contacts is proportional to q 2 /r 3 , while the pairwise interaction between 
monomers scales like l/(gp) (see eq [4]), giving [l3[] U ~ g/(pr 3 ) ~ l/(rb* 2 p). Therefore, for 

G(r,s) ~G^°)(r,a)(l-tO 

, <5( o,,, s) ( 1 _^), (27) 

which qualitatively explains eq[231 For q -j=^, we get: 

6 \ i 4tt 
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FIG. 7: For small internal end-to-end distances (r <C \/sb*), monomers of terminal blobs with g ~ r 2 /b* 2 
units interact. 

This regime is also limited by the condition qt; <C 1. Thus the low-g correction is positive and 
it increases with q, while the high-g correction 5G(q, s) is negative and is also increasing with q, 
implying an intermediate decline. This non-monotonic behavior of SG(q, s) translates in a non- 
monotonic dependence of 6F(q) for finite N (Figure [6]). 

To underline the origin of the "non-analytical" term 1/q, let us give some details of an easy 
calculation of the correction to the correlation function for the infinite chain. The top-left diagram 
in Figure [T] is the only diagram producing this term. With notations shown in Figure we can 
write the corresponding analytical expression: 

5F(q) = - J dku cff(k) exp(-k 2 (m - l)a 2 ) 

n~m,m-l , 0Q v 
t-n 

exp(— (q + k) 2 (n — m)a 2 ) exp(— k 2 (n — t)a 2 ). 

Here e - k2 ( m -0 a2 [ s the Gaussian (bare) correlation function of the chain between monomers I and 
m. In the interaction, we will leave only the fe 2 -dependent part. After summation (integration) 
over m — I and n — t we get: 

SF(q) = \- ^2 exp(-(q + k) 2 (m-n)a 2 )dcos(9)dk. (30) 

a n-m=0 •* 

Then, after integration over n — m we get: 

< 2 f 1 dcos(9) _ 1 1 
J-i Vl-cos 2 (0) ~ ~'pb*^]q\b* 
which is similar to eq [28j It should be emphasized that for a finite chain, the finite summation 
over m — I and n — t produces extra terms, and these "odd" terms exist only for qR„ < 1. 



«■<«> = -73b / ,. L, = -rknkt PD 




FIG. 8: Diagram producing the non-analytical term for an infinite chain. The mixing of the wave-vectors is 
pinned in the (m,n) segment. The summation is done on (n,t) and (m,l). 
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III. COMPUTATIONAL RESULTS 



The previous section contains non-trivial results due to generic physics which should apply to 
all polymer melts containing long and preferentially flexible chains. We have put these predictions 



261 ] of linear polymer melts having 



to a test by means of extensive lattice Monte Carlo simulations 
either a quenched and monodisperse or an annealed size-distribution. For the latter "equilibrium 
polymers" (EP) one expects (from standard linear aggregation theory) a Flory distribution if the 



27 



30|]. (Apart 



scission energy E is constant (assuming especially chain length independence) 
from this finite scission energy for EP all our systems are perfectly athermal. We set ksT = 1 and 
all length scales will be given below in units of the lattice constant.) We sketch first the algorithm 
used and the samples obtained and discuss then the intramolecular correlations as measured by 
computing the single chain form factor F(q). 



A. Algorithm and some technical details 



For both system classes we compare data obtained with the three dimensional bond-fluctuation 



model (BFM) 25J] where each monomer occupies the eig 



lattice. For details concerning the BFM see Refs. 



it sites of a unit cell of a simple cubic 



26, 



281 ] . For all configurations we use 
cubic periodic simulation boxes of linear size L = 256 containing 2 20 monomers. These large 
systems are needed to suppress finite size effects, especially for EP [28]. The monomer number 
relates to a number density p = 1/16 where half of the lattice sites are occupied (volume fraction 
0.5). It has been shown that this "melt densit y" is characterized by a small excluded volume 



correlation length £ of order of the monomer size 



25, 



261 ] . by a low (dimensionless) compressibility, 



S(q — * 0)/p = = 0.24, as may be seen from the bold line (/ = 1) indicated in the main panel 
of Figure [TOl below, and an effective statistical segment length, b* ~ 3.24 Q, [3]. 

The chain monomers are connected by (at most two saturated) bonds. Adjacent monomers are 
permanently linked together for monodisperse systems (if only local moves through phase space are 
considered). Equilibrium polymers on the other hand have a finite and constant scission energy E 
attributed to each bond (independent of density, chain length and position of the bond) which has 
to be paid whenever the bond between two monomers is broken. Standard Metropolis Monte Carlo 
is used to break and recombine the chains. Branching and formation of closed rings are explicitly 



forbidden and only linear chains are present. The delicate question of "detai 
microscopic reversibility, of the scission-recombination step is dealt with in Ref. 



ed balance", i.e. 



29] 
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28, 



29(|) that the 



The monodisperse systems have been equilibrated and sampled using a mix of local, slithering 
snake and double pivot moves extending our previous studies pj[ [l3| • This allowed us to generate 
ensembles containing about 10 3 independent configurations with chain length up to N = 8192 
monomers. More information on equilibration and sample characterization will be given elsewhere 
[49] ]. We have sampled EP systems with scission energies up to E = 15, the largest energy 
corresponding to a mean chain length (N) ~ 6000. Some relevant properties obtained for these 
EP systems are summarized in Table HI It has been verified (as shown in Refs. 
size-distribution obtained by our EP systems is indeed close to the Flory distribution studied in 
the analytical approach presented in the previous section. For EP only local hopping moves are 
needed in order to sample independent configurations since the breaking and recombination of 
chains reduce the relaxation times dramatically, at least for large scission-recombination attempt 
frequencies joll ]. In fact, all EP systems presented here have been sampled within four months 
while the sample of monodisperse N = 8192 configurations alone required about three years on the 
same XEON processor. EP are therefore very interesting from the computational point of view, 
allowing for an efficient test of theoretical predictions which also apply to monodisperse systems. 

B. Form factor 

It is well known that the intramolecular single chain form factor of monodisperse polymer chains 
may be computed using 

F N (q) = ^ (E • r *)) + (E sin (q • r *)) ) > ( 32 ) 

the average being taken over all chains and configurations available. The form factor obtained for 
the largest available monodisperse chain systems currently available is represented in the Figures [U 
El[7]and[9j It should be emphasized that the correct generalization of eq[32]to polydisperse systems 
compatible with eqs[6]and[T9]is the average with weight Npjy over F^(q). In practice, one computes 
simply the ensemble averaged sum over (X)i=i sin(q • rj) 2 + (J3i=i cos(q • rj) 2 contributions for 
each chain and divides by the total number of monomers. 

Figure [2] presents the (unsealed) form factors obtained for four different scission energies for 
our BFM EP model at p = 1/16. The three different g-regimes are indicated. Details of the 
size-distribution must matter most in the Guinier regime which probes the total coil size. Non- 
universal contributations to the form factor arise at large wavevectors (Bragg regime). Obviously, 
the larger E the wider the intermediate Kratky regime (see the dashed line indicating eq HOj) where 
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chain length, polydispersity and local physics should not contribute much to the deviations of the 
form factor from ideality. A very similar plot (not shown) has been obtained for monodisperse 
polymers. Not surprisingly, it demonstrates that the form factors of both system classes become 
indistinguishable for large wavevectors. 

The natural scaling attempt for the form factor of EP is presented in Figure [3] for a broad 
range of scission energies. We plot F(Q)/F(0) as a function of Q = qR g where both F(q — > 0) = 
{N 2 ^ / (N) and the (Z-averaged) gyration radius R g have been measured directly for each E. Note 
that the strong variation of F(0) and R g with E showing that the successful scaling collapse is 
significant. Obviously, this scaling does not hold in the Bragg regime (q ~ 1) where F(q) increases 
rapidly, as one expects. The bold line represents the ideal chain form factor, eq where the 
identification of the coefficients, F(0) — > 2/fj, and R g 2 — > 3a 2 / fj,, is suggested by the Guinier limit, 
eqEJ Hence, the perfect fit for q <C l/i? g is imposed, but the agreement remains nice even for much 
larger wavevectors. A careful inspection of the Figure reveals, however, that eq [7J overestimates 
systematically the data in the Kratky regime. (The corresponding plot for monodisperse chains is 
again very similar.) 

This can be seen more clearly in the Kratky representation given in Figure H] in linear coordi- 
nates. We present here the systems with the longest masses currently available for both monodis- 
perse (N = 4096 and N = 8192) and EP systems (E = 14 and E = 15). The non-monotonous 
behavior is in striking conflict with Flory's hypothesis. The difference between the ideal Gaussian 
behavior (thin line) and the data becomes up to 20%. For the large (average) chain masses given 
here all systems are identical for q > 0.1 (but obviously not on larger scales). It should be noted 
that qualitatively similar results have been reported — albeit for much shorter chains — for more 



than a decade in the literature 
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441 ] . The infinite chain prediction, eq[23](or equivalently eq 



gives a lower envelope for the data which fits reasonably — despite its simplicity — in the finite 
wavevector range 0.1 < q < 0.4. 

The form factor difference 5F(q) = F(q) — F^(q) is further investigated in the Figures [5] and [6] 
for equilibrium and monodisperse systems respectively. These plots highlight the deviations in the 
Guinier regime. In both cases the ideal chain form factor F^°\q) is computed assuming the same 
effective statistical segment length b* = a*^/6, i.e. the reference chain size is 

4°) = a *j\rV2, In 

the first case the reference F^°\q) is the ideal chain form factor for Flory-distributed chains, eq[JJ 
in the second the Debye function F^°(q) = NfniQ 2 ) with fo(x) = (exp(— x) — 1 + x)2/x 2 [3J. 

As suggested by eq [121 and eq 1201 respectively we plot 5F(q)/c vs. Q = qR^\ i.e. the axes 
have been chosen such that the data should scale for different (mean) chain length. We obtain 
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indeed a reasonable scaling considering that our chains are not large enough to suppress (for the Q- 
range represented) the deviations 5F(q) due to local physics. The scaling shows implicitly that the 
corrections with respect to the infinite chain limit decay as the inverse gyration radius, 1/c ~ 1/y/N, 
as predicted by eqs [18] and [21] (Both plots appear to improve systematically with increasing chain 
length and, clearly, high precision form factors for much larger chains must be considered in future 
studies to demonstrate the scaling numerically.) Also the functional agreement with theory is 
qualitatively satisfactory in both cases, for equilibrium polymers it is even quantitative for small 
wavevectors. For monodisperse chains we find numerically a much more pronounced hump in the 
Guinier as the one predicted by eq[20] (bold). This is very likely due to the chosen interaction 
potential eq H] for Flory-distributed chains which is not accurate enough for the description of the 
Guinier regime of monodisperse chains. 

It should be pointed out that the success of the representation of the non-Gaussian deviations 
chosen in the Figures [5] and [6] does depend strongly on the accurate estimation of the statistical 
segment length b* of the ideal reference chains. A variation of a few percents breaks the scaling and 
leads to qualitatively different curves. Since such a precision is normally not available (neither in 
simulation nor in experiment) it is interesting to find a more robust representation of the form factor 
deviations which does not rely on b* and allows to detect the theoretical key predictions for long 
chains (notably eq 1231) more readily. Such a representation is given in Figure [9] for monodisperse 
chains. (A virtually indistinguishable plot has been obtained for EP.) The reference chain size is 
set here by the measured radius of gyration R g (N) (replacing the above Rg 1 ^) which is used for 
rescaling the axis and, more importantly, to compute the Debye function F^°\q). The general 
scaling idea is motivated by Figure El the scaling of the vertical axis is suggested by eq [23] which 
predicts the difference of the inverse form factors to be proportional to N°q 3 . Without additional 
parameters (R g is known to high precision) we confirm the scaling of 

m(Q) s im - ms) plp " m 

as a function of Q = qR g with p* = N/R g 3 being the overlap density. Importantly, our simulations 
allow us to verify for Q ^> 5 the fundamentally novel Q 3 behavior of the master curve predicted 
by eq[23]and this over more than an order of magnitude! 

In this representation we do not find a change of sign for the form factor difference (5F(q) is 
always negative) and all regimes can be given on the same plot in logarithmic coordinates. In the 
Guinier regime we find now m(Q) oc Q 4 which is readily explained in terms of a standard expansion 
in Q 2 . (The first two terms in Q° and Q 2 must vanish by construction because of the definition of 
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FIG. 9: Scaling attempt of the non-Gaussian deviations of the form factor of monodisperse polymers in the 
melt in terms of the measured radius of gyration R g (instead of R^)- As suggested by eg I23l the difference 
1/F(q) — 1/F(°\q) of the measured and the ideal chain Debye form factor has been rescaled by the factor 
Np/p* and plotted as a function of Q = qR g . We obtain perfect data collapse for all chain lengths included. 
(Obviously, data points in the Bragg limit q w 1 do not scale.) The difference —SF(q) is positive in all 
regimes and no change of sign occurs in this representation. Note that the power law slope, m(Q) = Q 3 /32, 
predicted by eq I23L can be seen over more then one order of magnitude. In the Guinier regime, the difference 
increases more rapidly, m(Q) oc Q (dashed line), as one expects from a standard analytic expansion in Q 2 . 

radius of gyration, eq[8] [3].) Finally, we stress that the scaling of Figure [9] is not fundamentally 
different from the one attempted in Figures [5] and [6J Noting F(q)F^ (q) « (N fo(Q)) 2 , it is 
equivalent to -8F(Q)/c ra m(Q)f D (Q) 2 with c w N 2 /pR g 3 . (Compared with eqs QJ and EQl 
has been replaced by the measured R g .) This scaling has been verified to hold (not shown) but we 
do not recommend it, since it does not yield simple power law regimes. 



22 




FIG. 10: The response function S(q, /) of a fraction / of marked monodisperse chains of chain length N = 
512. The remaining 1 — / chains are considered to be not "visible" for the scattering. (The simulation box of 
linear size L = 256 contains 2048 chains in total.) The main figure (a) presents S(q, /) directly as computed 
from S(q, f)/pf — ^ ((^2i cos(q • r^)) 2 + sm (q ■ rj) 2 ^} where the sums run over all n — L 3 pf marked 
monomers and the wavevectors are commensurate with the cubic box. Also included is the form factor F(q) 
(dashed line) from eg [521 which corresponds to the S(q, f)/pf F(q) limit (but compares already perfectly 
with the / = 0.01 data set). The so-called "total structure factor" S(q,f — 1) (bold line at bottom) is the 
Fourier transformed monomer pair-correlation function of all monomers. For scales larger than the monomer 
size it is essentially constant and yields an accurate determination of the (dimensionless) compressibility 
S(q)/p — 1/vp « 0.24. Since the compressibility is finite, although small, this value has to be properly 
taken into account if one wants to estimate the form factor F(q) from the experimentally available S(q, /). 
The inset (b) presents a Kratky representation of S(q, /) as suggested by an appropriate generalization of 
eqffl S(q, /) - } 2 S(q, 1) = pf(l - f)F(q); the reduced quantity q 2 [S(q, f) - f 2 S(q, 1)] / [pf(l - /)] Q, Q 
is plotted against the wavevector. In addition, q 2 F(q) obtained using eq[32]is indicated by the dashed line 
in the inset. All data sets collapse perfectly which provides a striking confirmation of eq[T] Also indicated 
is the infinite chain asymptote eq[22] (bold line). 
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IV. CONCLUSION 

We have shown in this paper that even for infinitely long and flexible polymer chains no Kratky 
plateau should be expected in the form factor measured from a dense solution or melt (see Figured]). 
We rather predict a non-monotonic correction to the ideal chain scattering crossing from positive in 
the Guinier regime to negative in the Kratky regime (Figures [5] and [6]) . The former regime merely 
depends on the radius of gyration and the correction corresponds to some deswelling of the coil. In 
the latter regime the form factor ultimately matches that of an infinite chain = ^2 — 1~ jil^r 
for q£ <C 1 (Figure M)- 

The t/ 3 -correction depends neither on the interaction nor on the statistical segment, it must 
hence be generally valid, even in the critical semidilute regime. We checked explicitly that the 
one-loop correction obtained by Schafer in the strongly fluctuating semidilute regime by numerical 
integration of renormalization group equations takes the same form with the amplitude 0.03124 .. . 
within 0.03% of our 1/32. 

It is to be noted that the above correction for infinite chains is not an analytic function of q 2 
as one would naively anticipate. For finite chains the correction remains a function of even powers 
of q. The intriguing g 3 -correction for infinite chains formally arises from dilation invariance of the 
diagrams. Established theoretical methods 
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52l | may implicitly assume analytical properties 



of scattering functions and non-analytical terms discussed in our paper could be easily overlooked. 

These theoretical results are nicely confirmed by our Monte Carlo simulations of long flexible 
polymers. The agreement is particulary good for equilibrium polymers (Figure [5]) and satisfactory 
for all systems with large (mean) chain length N 1000 (Figure ED . It should be emphasized that 
all fits presented in this paper are parameter free since the only model dependent parameter b* has 
been independently obtained from the internal distances of chain segments [48]. Since a sufficiently 
accurate value of b* may not be available in general, our simulation suggest as a simple and robust 
way to detect (also experimentally) the universal g 3 -correction the scaling representation of the 
(inverse) form factor difference in terms of the measured radius of gyration given in Figure [9J 
We expect that data for any polymer sample — containing long and flexible linear chains with 
moderate polydispersity — should collapse with good accuracy on the same master curve. Strong 
polydispersity (such as one finds in EP) should merely change its behavior in a small regime around 
Q = qR g K, 1. 

Measuring the form factor is a well accepted method to determine the statistical segment length. 
We already mentioned in the Introduction that there is no clear evidence of a true Kratky plateau 
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from experiments and further showed that a plateau is actually not to be expected from the theory 
on general grounds. We are lacking an operational definition of the statistical segment length, even 
for very long flexible "thin" chains. One way out would be in principle to fit a large g-range, from 
the Guinier regime — as far it can be cleanly measured on a sample with controlled polydispersity 
- to monomer scale, with the corrected formula (q) + 5F(q). However, if the size-distribution 
is not known precisely (as it will be normally the case) we recommend to determine b* instead by 
means of the infinite chain asymptote, eq 1221 as can be seen Figure HI 

At this point one may wonder whether eqQ3 (the precise form of this equation being given in 
the caption of Figure [TO]) routinely used to rationalize the scattering of a mixture of deuterated 
and hydrogenated chains is accurate enough to extract the form factor, including the corrections. 
From a theoretical point of view, for "ideal" labeling of the chains, which does not introduce 
additional interactions between labeled and unlabeled chains, there is no question that this can be 
done. Practically however, there is a danger that experimental noise in subtracted terms in eq[T] 
will mask corrections discussed in our paper. The strongest support comes here from numerical 
results presented in Figure [TOl We have computed the response function S(q, f) for a melt of 
monodisperse chains for chain length N = 512 and different fractions / of labeled chains. The 
main panel (a) gives S(q, /)/ fp and the form factor F(q) as a function of the wavevector. The 
inset (b) presents a Kratky representation of the rescaled structure factor: For a surprisingly large 
range of / the data scales if the standard experimental procedure is followed and the scattering 
of the background density fluctuations f 2 S(q,f = 1) has been properly substracted [j], 43]. The 
rescaled response function is identical to F(q) and shows precisely the non-monotonic behavior, 
eqEH and the asymptotic infinite chain limit, eq[22] (bold line in inset), predicted by our theory. 
This confirms that eq Q] allows indeed to extract the correct form factor and should encourage 
experimentalists to revisit this old, but rather pivotal question of polymer science 3]. 
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143 


15 


6183 
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TABLE I: Various properties of linear equilibrium polymers obtained with the BFM algorithm at number 
density p = 0.5/8 (volume fraction 0.5): imposed scission energy E, the mean chain length (A 7 ), the ratio 
F(0) — (N 2 ) /(N) comparing the first and the second moment of the number distribution, the mean-squared 
bond length I 2 = (l 2 ), the z-averaged end-to-end distance R e and the radius of gyration R g . For all scission 
energies we have used periodic simulation boxes of linear size L = 256 containing n mon = 2 20 monomers. 
Averages are performed over all chains and 1000 configurations. 
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as indicated by eqs [T5] and [TH] for the polydisperse case. Due to chain end effects it turns out that this 





28 



requires much larger (mean) chain lengths than the recommended method above. 
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Since the EP chains break and recombine permanently the relaxation time of the system is not set by 
the typical EP radius of gyration but rather by the size of a small chain segment which has just about 
the time to diffuse over its radius before it breaks or recombines 27]. The high frequency for the scission- 
recombination attempts used in our simulations ensures that the effective recombination time is small 
and the dynamics is, hence, always of Rouse-type. It should be emphasized that due to the permanent 
recombination events a data structure based on a topologically ordered intra chain interactions is not 
appropriate and straight-forward pointer lists between connected monomers are required [28j. The 
attempt frequency should not by taken too large to avoid useless immediate recombination of the same 
monomers and some time must be given for the monomers to diffuse over a couple of monomer diameters 

n 

between scission- recombination attempts [29] . 
Fuchs, M., Z. Phys. B 1997, 521-530. 

Schulz, M.; Frisch, H.L.; Reineker, P., A^ew Journal of Physics 2004, 6, 77. 

We have computed the response function S(q, f = 1) for monodispcrse polymers over a large range of 
densities and for a nouvel version of the BFM with finite overlap energies. This has been done to verify 
the recent prediction of fluctuation induced long-range repulsions between solid objects in polymer 



media [T5L il 6t| . This approach suggests a systematic violation of the RPA eq [5] proportional to q 3 which 
we have put to a numerical test. Our findings — complicated by the fact that trivial monomer- monomer 
correlations of Percus-Yevick type 
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must be correctly taken into account |24j — will be presented 
elsewhere. Please note that these corrections to eq [2] do correspond to higher order deviations which 
can be shown to be negligible for the questions addressed in this paper. 
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